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ABSTRACT 

In  publicly  available  generic  digital  simulations  of  semi-active  Radio  Frequency 
(RF)  seekers,  the  combined  effect  of  multipath,  sea-clutter  and  other  non¬ 
linear  noise  sources  are  mainly  studied  in  individual  models,  while  the  effect 
in  the  seekers  data  processing  has  not  been  systematically  investigated.  This 
report  details  these  effects  and  how  they  have  been  combined  to  model  the 
signal  processing  blocks  inside  a  semi-active  seeker.  A  basic  doppler  detection 
scheme  is  then  simulated  to  investigate  the  consequence  of  these  effects  on 
target  detection. 
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Signal  Processing  in  a  Semi-active  Seeker 


EXECUTIVE  SUMMARY 

In  a  generic  semi-active  RF  seeker  in  medium  or  long  range  missiles,  the  combined 
effect  of  multipath,  sea-clutter  and  other  non-linear  noise  sources  may  cause  detrimental 
effects  on  target  detection.  This  report  describes  models  for  each  of  the  effects  mentioned 
and  how  they  can  be  combined  into  a  single  model.  A  simulation  of  the  combined  effects 
is  then  used  in  a  simple  detection  scheme,  so  that  any  change  in  target  detection  can  be 
assessed. 
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1  Introduction 


1.1  Project  Aim 


In  a  generic  semi-active  RF  seeker  in  medium  or  long  range  missiles,  the  signals  at 
the  front  and  rear  receivers  are  mixed  to  shift  the  carrier  frequency  to  an  Intermediate 
Frequency  (IF)  where  the  target’s  doppler  shift  can  be  detected.  When  the  missile  is  used 
against  low  altitude  targets  in  the  presence  of  multipath  in  the  rear  receiver  and  sea-clutter 
and  multipath  in  the  front  receiver,  the  mixing  of  these  two  signals  will  significantly  affect 
the  seeker’s  ability  to  detect  a  target  correctly. 


In  publicly  available  generic  digital  simulations  of  semi-active  RF  seekers,  there  has 
been  substantial  research  into  both  multipath,  [l]-[3]  and  sea-clutter,  [4]-[5].  The  aim  of 
this  project  is  to  combine  these  effects  into  a  single  model,  simulate  them  and  investigate 
the  mixing  of  the  front  and  rear  signals.  A  target  detection  scheme  is  then  used  to  deter¬ 
mine  the  effect  on  target  detection. 


The  report  starts  by  describing  a  classical  model  for  a  Continuous  Wave  (CW)  seeker 
and  a  simplified  version  used  for  simulation.  The  second  section  describes  three  effects: 
multipath,  sea-clutter  and  receiver  noise,  which  are  used  to  increase  the  realism  in  the 
model.  The  final  section  focusses  on  doppler  detection  and  outlines  an  optimal  detection 
scheme  as  well  as  a  classical  non-optimal  detector  described  in  Skolnik,  [6] .  This  is  exten¬ 
sively  modelled  and  Receiver  Operating  Characteristics  (ROC)  are  compared  for  both  an 
ideal  simulation  and  a  realistic  one  using  the  three  effects. 


1.2  Problem  Description 

A  typical  encounter  with  a  missile  containing  a  semi-active  RF  seeker  against  sea- 
skimming  targets  is  shown  in  Figure  1.  The  seeker  system  has  two  receivers,  one  at  the 
front  and  the  other  at  the  rear.  As  a  semi-active  system,  it  contains  no  active  radar  and 
is  completely  dependent  on  the  launch  platform  to  illuminate  the  target  and  provide  rel¬ 
evant  information  about  the  frequency  of  the  illuminating  beam  and  its  modulation,  if  any. 


As  the  target  and  the  missile  approach  each  other,  two  doppler  shifts  arise.  The  first 
is  from  the  launch  platform  to  the  missile  and  the  second  is  from  the  launch  platform 
reflected  off  the  target  and  incident  on  the  seeker.  This  differential  doppler  allows  the 
seeker  to  identify  and  lock  onto  the  target  by  carefully  combining  and  filtering  the  two 
incoming  signals.  To  calculate  the  position  of  the  target,  a  monopulse  receiver  at  the 
front  of  the  seeker  measures  the  angle  of  arrival  (AOA)  of  the  target  signal  in  azimuth  and 
elevation.  This  information  will  be  noisy  however,  due  to  multipath  reflections,  sea-clutter 
and  non-linearities  from  the  hardware  in  the  seeker. 
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Reflections  off  sea  surface 


Figure  1:  Problem  scenario  with  multipath  and  sea-clutter  reflections 

2  The  Semi- Active  Seeker 


2.1  Classic  Model 

The  block  diagram  in  Figure  2,  [6],  is  representative  of  the  earliest  systems  developed 
in  the  late  1940’s  and  early  1950’s.  It  is  the  simplest  representation  of  a  semi-active  CW 
seeker  and  consists  of  a  rear  receiver,  a  front  receiver,  a  signal  processor  (speedgate),  and  a 
tracking  loop  to  control  the  gimballed  front  antenna.  The  missile  also  contains  an  autopi¬ 
lot  to  guide  it  and  stabilize  the  airframe,  a  fuze  to  detonate  the  warhead  at  the  optimum 
time,  and  a  source  of  electrical  and  (in  most  missiles)  hydraulic  power. 
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The  purpose  of  the  rear  reference  receiver  is  to  provide  a  coherent  reference  for  the 
detection  of  the  front  or  target  signal.  The  rear  signal,  after  conversion  to  IF,  closes  the 
Automatic  Frequency  Control  (AFC)  loop  around  the  microwave  Local  Oscillator  (LO) 
and  acts  as  the  reference  for  the  IF  coherent  detector.  The  target  signal,  received  at  the 
front  antenna,  is  heterodyned  to  IF  and  amplified  in  a  relatively  wideband  amplifier.  It  is 
then  converted  to  the  baseband  by  mixing  with  the  rear  signal  in  the  balanced  mixer. 


The  doppler  signal,  now  at  baseband,  is  amplified  in  the  video  (doppler)  amplifier, 
which  has  a  bandwidth  equal  to  the  total  range  of  possible  doppler  frequencies.  It  is  then 
mixed  with  the  speedgate  LO,  which  is  controlled  by  an  AFC  loop  to  keep  the  desired 
signal  centered  in  the  narrow  speedgate  (sometimes  called  the  velocity  gate  or  doppler 
tracker).  Target  detection  is  covered  more  thoroughly  in  section  4. 


2.2  Simplified  Model 

The  semi-active  seeker  model  detailed  in  the  previous  section  is  now  simplified  so  the 
effect  of  different  noise  sources  in  the  doppler  output  can  be  studied. 

•  The  AFC  has  been  removed  as  both  front  and  rear  signals  are  created  from  ideal 
frequency  sources. 

•  The  gimbal  servo  and  antenna  control  has  been  removed,  as  the  front  receiver  is 
assumed  to  point  to  the  target  and  the  rear  receiver  is  fixed  at  the  exit  direction 
with  a  broad  view  angle. 

•  The  first  mixing  stage  including  the  LO  has  also  been  removed,  as  simulating  an 
RF  signal  in  the  Gigahertz  band  is  not  feasible.  Also,  the  effect  of  amplifying  the 
IF  signals  in  the  Front  and  Rear  IF  blocks  may  not  have  a  big  effect  on  the  final 
doppler  measurement.  A  low  pass  filter  has  been  included  in  its  place  and  is  placed 
after  the  mixer  to  remove  unwanted  high  frequency  signals. 

•  The  section  after  the  balanced  mixer  is  addressed  separately  in  section  4. 

The  resulting  model  is  shown  in  Figure  3. 

The  mixing  of  the  front  and  rear  signals,  Sr  and  Sp  is  either  a  multiplication  in  the 
time  domain,  or  a  convolution  in  the  frequency  domain. 

In  the  time  domain: 

Soit)  =  SR{t)SF{t)  (1) 

In  the  frequency  domain: 

SD{f)  =  SR{f)^SF{f)  (2) 


The  low  pass  filter  is  then  used  to  remove  unwanted  high  frequency  components  (see 
section  2.3). 
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Spectral  Measwement 


Figure  3:  Simplified  semi-active  seeker  signal  processing  block  diagram 

2.3  Doppler 

The  doppler  effect  is  a  shift  in  the  frequency  of  a  wave  radiated,  reflected  or  received  by 
an  object  in  motion.  In  the  scenario  for  this  project,  a  stationary  illuminator  is  tracking  a 
moving  target  and  the  radio  waves  either  compress  or  stretch  when  they  come  in  contact 
with  the  target.  The  doppler  frequency  equals  the  rate  of  change  of  the  distance  to  the 
target  d,  divided  by  the  wavelength. 

fd  =  (3) 

The  minus  sign  accounts  for  the  fact  that,  if  d  is  negative  (closing  target),  the  doppler 
frequency  is  positive. 


The  seeker  makes  use  of  the  differential  doppler  frequency  to  distinguish  small  but  fast 
moving  airborne  targets  from  each  other  and  from  any  background  terrain  and  clutter. 
Early  doppler  homers  used  pure  CW  illumination  of  the  target,  but  recently  the  tendency 
has  been  to  employ  some  type  of  modulation  such  as  pulse  or  CW  frequency  modulation. 
A  system  with  a  modulation  scheme  is  able  to  determine  both  range  and  range  rate 
simultaneously,  but  for  simplicity  the  simulation  in  this  project  is  a  pure  CW  system.  The 
geometry  for  an  ideal  doppler  shift  of  the  front  and  rear  receivers  is  shown  in  Figure  4, 

[7]. 


Let  the  illuminator  (I)  carrier  frequency  be  /c)  the  missile  and  target  velocities,  Um  and 
ut  respectively.  The  target  (T)  receives  a  frequency  /c  +  {^)cosfi  due  to  a  movement  of 
T  with  respect  to  I.  This  then  re-radiates  towards  the  missile  (M)  with  a  fiirther  doppler 
shift  of  (^)  cos  a  due  to  M  travelling  towards  T.  The  missile  is  also  affected  by  that  motion 
and  has  a  doppler  shift  of  {^)cos6.  Hence  the  front  receiver  on  the  seeker  receives  a 
signal  with  frequency 


fc  + 


Um  cos  S  +  Mt(cos  a  +  cos  0) 
A 


(4) 


The  mixer  inside  the  seeker  extracts  the  target’s  doppler,  /d,  by  comparing  this  frequency 
with  that  of  the  rear  reference,  which  is 


fc- 


Um  cos  J 

A 


(5) 
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Figure  4-'  Semi-active  homing  target  geometry 


due  to  a  movement  of  M  away  from  I.  Now,  define 

cos  5  +  Wt  (cos  O' -f- cos /?) 

A 

Um  COS  7 

A 

and  when  two  CW  signals  at  these  frequencies  are  multiplied  in  the  time  domain  as  in 
equation  1,  the  following  occurs; 


fp 

fp 


SD{t)  =  Ap sm{27r(f c fF)t)ARsm{2-K{fc  +  fR)t) 

=  — ^ —  [cos(27r(/ir  -  fR)t)  -  cos(27r(2/c  +  fp  +  fnW  (6) 

where  each  signal  has  zero  phase.  After  low-pass  filtering  the  result,  and  sampling  high 
enough  to  see  only  the  bottom  part  of  the  spectrum  where  the  final  doppler  will  be,  the 
final  amplitude  of  the  doppler  signal,  Sp  is: 


Ad 


ApAr 

2 


with  corresponding  doppler  frequency: 


(7) 


Id  =  fp  -  Jr 


Um  (cos  5)  -h  Wf  (cos  a  -1-  cos  pi)  ^  Um  cos  7 

A  ^  3r“ 

Um{cos  5  -I-  COS  7)  4-  Wf  (cos  a  -h  cos  P) 

A  ■ 


(8) 


The  maximum  value  of  fp  occurs  when  I,  M  and  T  are  collinear  and  M  and  T  are  flying 
towards  each  other,  giving 


2(^771  T  'ut) 

T  ■ 


(9) 
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3  Simulated  Effects 

Ekjuation  8  represents  ideal  line  spectra.  However,  in  real  life  the  spectra  spreads  in 
frequency,  making  computation  more  complicated.  Effects  such  as  main  lobe  sea-clutter 
which  will  cause  an  extra  peak,  multipath  reflections  which  spread  the  front  and  rear  sig¬ 
nals  and  receiver  noise  will  alter  the  ideal  spectrum.  To  explore  these  effects,  simulation 
is  required. 

The  simulation  for  this  section  was  written  in  MATLAB,  using  Simulink  for  the  sim¬ 
ulation  and  signal  processing.  The  CW  or  carrier  frequency  is  chosen  as  lOGHz  which 
corresponds  to  an  X-band  radar.  However,  the  concerned  positive  and  negative  doppler 
effects  will  be  around  lOGHz.  Therefore,  the  CW  frequency  was  used  as  a  reference  and 
a  sampling  rate  cover  range  of  4MHz  was  used  to  give  a  spectra  range  of  10GHz±2MHz. 


The  simulation  is  done  entirely  in  the  frequency  domain  to  simplify  the  addition  of  the 
external  effects.  Consequently,  the  mixing  is  now  a  convolution  instead  of  a  multiplication 
in  the  time  domain.  The  effect  of  convolving  two  signals  is  a  lengthening  of  the  vector, 
a  final  amplitude  dependent  on  the  length  of  the  input  and  triangulation.  For  example, 
if  two  uniform  signals  Xi  and  X2  have  lengths  JV  and  M  respectively,  then  Xi  ®  X2  has 
length  M  -|-  iV  —  1  and  the  spectrum  is  shown  in  Figure  5. 


X,  X3 


Figure  5:  Convolution  example 


The  CW  signal  at  the  illuminator  is  modelled  as  a  Gaussian  shape  with  a  small  fre¬ 
quency  spread  and  a  mean  frequency,  pcw  equal  to  lOGHz. 


Scwif)  = 


:exp 


27r<r: 


CW 


if  -  l^cwf 


2o'qw 


The  front  and  rear  CW  signals  are  then  created  in  the  simulation  by  adding  a  frirther 
doppler  shift  specified  by  the  parameters  in  the  scenario.  The  frequency  offsets  for  the 
front  and  rear  signals  are  defined  as  pi  and  p2- 
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To  increase  realism,  three  noise  sources  have  been  included  before  mixing,  (1)  multipath 
in  the  front  receiver,  (2)  main-lobe  sea-clutter  in  the  front  receiver  and  (3)  receiver  noise 
in  both  the  front  and  rear  receivers.  The  simulation  runs  for  a  single  time  period,  so  the 
missile  and  the  target  are  at  a  constant  velocity  and  do  not  change  their  position.  Two 
test  scenarios  were  created  for  the  analysis  and  have  the  following  parameters: 


Table  1:  Scenario  parameters 


Scenario 

no. 

a(“) 

13  n 

7(°) 

Missile  velocity, 
Um  (m/s) 

Target  velocity, 
ut  (m/s) 

Doppler  - 
front  (kHz) 

Doppler  - 
rear  (kHz) 

1 

0 

0 

0 

0 

900 

300 

50 

-30 

2 

60 

90 

60 

60 

900 

300 

20 

-15 

The  first  set  of  parameters  was  chosen  as  an  ideal  case  where  the  target,  missile  and 
illuminator  are  in  a  straight  line.  The  missile  velocity  is  set  to  three  times  the  target 
velocity  and  the  target  is  assumed  to  be  incoming.  The  second  set  of  parameters  are  a 
little  more  realistic  and  have  the  target  flying  straight  up  trying  to  avoid  the  incoming 
missile. 


3.1  Ideal  Simulation 

The  ideal  simulation  results,  without  any  extra  effects,  are  shown  in  Figures  6  and  7. 
They  are  based  on  equation  11,  which  is  represented  in  the  frequency  domain. 

SD  =  S)a(8)SF  (11) 

For  both  scenarios,  the  results  clearly  show  the  target’s  doppler  at  80kHz  and  35kHz 
respectively.  A  verification  for  this  simulation  is  given  in  Appendix  A. 


ideal_case1 


Frequency  (kHz) 

Figure  6:  Ideal  scenario  1 
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100  120 
Frequency  (kHz) 


Figure  8:  Rear  multipath  components 


A  duct  is  a  layer  of  air  with  different  properties  than  its  surroundings  and  acts  as  a 
waveguide  to  trap  electromagnetic  energy.  These  refractive  gradients  in  the  atmosphere 
can  cause  over-the-horizon  fields  to  be  tens  of  decibels  higher  than  expected.  More  detail 


Figure  1:  Ideal  scenario  2 


3,2  Multipath  Reflections 


Multipath  effects  come  about  by  reflections  on  the  surface  of  the  sea  and  by  ducting, 
[2].  Multipath  reflections  from  the  sea  can  be  either  specular  or  diffuse.  The  specular 
component  is  coherent  with  respect  to  the  direct  signal,  and  is  well  defined  in  terms  of 
amplitude,  phase  and  incident  direction.  The  diffuse  component  has  a  random  nature,  and 
arises  from  scattering  sources  from  many  directions.  These  paths  are  shown  in  Figure  8,  [3]. 


arUlrflNATOR 
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about  ducting  can  be  found  in  [2]  and  [6].  For  this  project,  only  the  effect  of  optical 
interference  has  been  included. 


Due  to  the  small  grazing  angle  in  the  scenario  for  this  project,  multipath  reflections  are 
likely  to  be  strong,  even  with  high  sea  states.  According  to  plane  earth  reflection  theory, 
the  one-way  voltage  reflection  factor  is  2sin(^2j^),  where  hi  and  h2  are  the  radar  and 
target  heights  respectively  and  R  is  the  range.  Practical  values  and  a  more  comprehensive 
explanation  are  given  by  James  in  [7]. 


In  the  semi-active  scenario,  there  are  three  ways  for  multipath  forward  scattering:  (1) 
the  transmit  path  from  illuminator  to  missile  rear  receiver,  (2)  the  transmit  path  from 
the  illuminator  to  target,  and  (3)  the  receive  path  from  the  target  to  missile  in  the  front 
receiver  (see  Figure  9,  [3]).  In  this  simulation,  we  consider  only  multipath  corruption  of 
the  rear  signal. 


Figure  9:  Front  multipath  components 

It  has  been  shown,  [1]  that  the  effect  of  multipath  reflections  distorts  the  final  doppler 
signal  by  widening  the  spectrum.  For  the  simulation  in  this  project,  the  rear  signal,  Sr 
will  change  from  a  line  spectrum  to  a  broad  spectrum,  where  the  spread  in  rear  signal 
accounts  for  both  the  specular  and  diffuse  paths. 

Multipath  in  the  front  receiver  was  studied  as  a  special  case  of  the  multiple  target 
problem,  where  resolving  two  closely  spaced  targets  is  very  difficult.  In  tracking  radar, 
the  presence  of  unresolved  targets  will  degrade  the  quality  of  the  data  obtained.  In  a 
guidance  system,  the  approach  used  to  overcome  this  is  to  minimize  the  reflected  (image) 
signal  and  also  prevent  the  ‘noise’  in  the  pitch  (elevation)  plane  guidance  channel  from 
causing  the  missile  to  impact  the  surface  during  low-altitude  intercepts.  In  the  case  of 
the  sea-skimming  anti-ship  missile,  this  is  done  by  implementing  a  radar  altimeter  for 
guidance  in  the  elevation  plane,  and  by  only  using  the  target  angle  data  in  azimuth. 

3.2.1  Simulation  Results 

There  are  many  models  that  can  be  used  to  simulate  the  multipath  signal.  The  mul¬ 
tipath  signal  for  this  project  is  created  by  the  original  CW  signal  with  no  doppler  shift. 
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Sew,  convolved  with  a  multipath  ftinction;  M.  M  is  a  Gaussian  centered  at  the  rear 
doppler  and  the  result  of  this  convolution  is  another  Gaussian  centered  with  a  frequency- 
shift  /i2,  the  same  as  the  rear  signal.  The  resultant  is  modelled  with  amplitude  Am  and 
a  very  broad  user  specified  frequency  spread  trl^.  In  reality,  the  spectrum  would  skew  to 
the  left  as  the  doppler  shift  in  the  diffuse  path  is  smaller. 


Suif)  =  SGw{f)®M{f) 


Am 


exp 


(/  - 


(12) 


The  rear  signal  can  then  be  repr^ented  as  the  sum  of  direct  and  diffuse  components: 


Sk  =  SR  +  SM  (13) 

The  spread  rear  signal  is  then  further  spread  by  the  convolution  and  the  equation  for  the 
doppler  output  signal  can  now  be  updated: 


Sd  =  (Sr  +  Sm)  ®  Sp 

Figures  10  and  11  demonstrate  the  multipath  reflections  for  each  scenario. 


(14) 


multipafri_Gas©1 


Figure  10:  Scenario  1  with  multipath  reflections 


mulUpafri_case2 


Figure  11:  Scenario  2  with  multipath  reflections 

In  each  of  the  above  figures,  it  is  clear  that  the  original  signal  is  greatly  spread  in 
frequency.  This  amount  of  spread  was  chosen  to  give  a  better  comparison  to  real  life. 
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3.3  Sea  Clutter 

As  the  sea-skimming  target  is  being  illuminated,  there  is  a  large  amount  of  background 
clutter.  This  sea  clutter  is  then  seen  by  the  front  receiver  on  the  missile  which  is  looking 
down  at  the  target  (see  Figure  12). 


Figure  12:  Front  and  side-lobe  clutter  for  a  sea-skimming  missile 


The  return  from  the  sea  is  made  up  of  two  main  parts,  main-lobe  and  side-lobe  return 
and  is  typically  different  for  different  sea  states  and  atmospheric  conditions.  In  general, 
the  power  of  the  return  from  a  small  patch  of  ground,  at  a  given  transmitter  frequency  is 
given  by 


where 


oc 


PavgG  CPAg 


(15) 


avff 


G 


cr  = 


■^9 

R 


average  transmitted  power 

gain  of  radar  antenna  in  the  direction  of  the  patch  (G^  =  two  way  gain) 
factor  called  the  incremental  backscattering  coefficient 
resolvable  area  of  ground  (ground  patch) 
range  of  ground  patch 


The  backscattering  coefficient,  a  is  the  radar  cross  section  of  a  small  increment  of  ground 
area,  A  A.  This  parameter  is  explained  further  in  [8]. 


3.3.1  Main-lobe  return 

Main-lobe  return  or  Main-lobe  Clutter  (MLC)  is  produced  whenever  the  main-lobe 
intercepts  the  sea.  Because  the  area  intercepted  by  the  main-lobe  can  be  extensive  and 
the  gain  of  the  main-lobe  is  high,  main-lobe  return  is  generally  quite  strong. 

The  spectral  characteristics  of  main-lobe  return  are  best  understood  by  visualizing 
the  ground  area  illuminated  by  the  main-lobe  as  consisting  of  a  large  number,  of  small 
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individual  patches.  The  doppler  frequency  of  each  patch,  fd  is  proportional  to  the  cosine 
of  the  angle  a,  between  the  missile  velocity  and  the  Line  of  Sight  (LOS)  to  the  patch. 

(16) 

where 


fd  — 


Uffi  cos  Q( 


Um  =  velocity  of  missile 

a  =  angle  between  the  missile  and  LOS  to  the  ground  patch 
A  =  wavelength 

The  angle  a  is  not  the  same  for  every  part  of  the  patch.  As  a  result,  the  collective  return 
occupies  a  band  of  frequencies. 

For  the  case  when  the  antenna  is  looking  straight  ahead,  the  doppler  frequency  of  the 
return  from  the  patches  near  the  center  of  the  illuminated  area  very  nearly  equals  its 
maximum  possible  value: 

(17) 


fd  = 

Mfnax 


llri 


3.3.2  Side-lobe  return 

The  clutter  signal  received  from  the  antenna  side-lobes,  also  called  Side  Lobe  Clutter 
(SLC)  is  always  undesirable.  It  tends  to  be  less  concentrated  (less  power  per  unit  of  doppler 
frequency)  than  the  main-lobe  clutter,  but  covers  a  much  wider  band  of  frequencies.  Side- 
lobes  extend  in  virtually  all  directions,  therefore  regardless  of  the  antenna  look  angle, 
there  are  always  side-lobes  pointing  ahead,  behind  and  every  angle  in  between.  The 
band  of  frequencies  covered  by  the  side-lobe  clutter  extends  from  a  positive  frequency, 
corresponding  to  the  radar  velocity  {fa  =  ^),  to  an  equally  negative  frequency.  The 
extent  to  which  side-lobe  clutter  is  a  problem,  depends  on  many  things: 

•  Frequency  resolution  provided  by  the  radar 

•  Range  resolution  provided  by  the  radar 

•  Gain  of  the  side-lobes 

•  Altitude  of  the  radar 

•  Backscattering  coefficient  and  angle  of  incidence 

•  Man-made  objects  in  the  terrain 

3.3.3  Simulation  Results 
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In  the  simulation,  only  MLC  has  been  included,  as  the  effect  of  SLC  on  doppler  is 
negligible  compared  to  a  peak  that  is  near  or  overlapping  a  target  signal.  The  sea-clutter 
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signal  is  created  in  the  frequency  domain  with  a  Gaussian  shape,  user  defined  amplitude 
Ac,  broad  frequency  spread  cr^  and  frequency  offset  fd- 


Seif)  = 


if -fcl)^' 


This  distribution  is  not  an  exact  representation  of  real  life,  as  there  are  many  parameters 
involved  which  shape  the  MLC.  It  was  chosen  as  a  simple  model  to  use  in  this  simulation 
and  could  be  improved  in  future  work.  In  the  simulation,  the  sea-clutter  signal  is  added 
to  the  front  signal,  which  now  becomes: 


Sp  =  Sf  +  Sc  (19) 

where  Sp  is  the  original  front  signal.  The  equation  for  the  doppler  output  signal  can  now 
be  updated; 

Sd  =  Sr  (g)  (Sp  -t-  Sc)  (20) 

where  the  effect  of  the  convolution  is  to  further  spread  the  sea-clutter  signal.  For  the 
scenarios  in  this  project,  the  desired  frequency  offset,  fd  is  set  to  the  corresponding  doppler 
between  the  missile  and  the  ground  offset  by  20°  to  represent  the  peak  sea-clutter  signal 
not  being  directly  at  the  line  of  sight. 


fd 


Um  cos  (20°) 
A 


(21) 


For  both  scenarios,  the  missile’s  velocity  Um  is  300m/s,  and  equation  21  gives  a  fre¬ 
quency  of  28.21kHz.  The  mixing  then  causes  the  center  frequency  to  shift  to  the  right  by 
subtracting  the  rear  doppler  component.  Hence,  the  expected  frequency  for  scenario  1  is 
28.21kHz  —  (—30kHz)  =  58.21kHz  and  28.21kHz  —  (—15kHz)  =  43.21kHz  for  scenario  2. 
This  calculation  is  detailed  in  Appendix  B. 


Figures  13  and  14  show  the  results  from  the  Simulink  model.  The  amplitude  of  the 
sea-clutter  signal  is  set  to  be  slightly  lower  than  the  target  signal,  and  the  spread  is  large 
to  give  a  better  comparison  to  real  life.  In  each  case,  it  is  clear  that  the  added  signal  is  at 
the  correct  frequency  and  the  characteristic  of  the  sea-clutter  signals  are  correct. 
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sea_clutter_case1 


Figure  13:  Scenario  1  with  sea-clutter 


sea_clutt©f_case2 


Figure  14:  Scenario  2  with  sea-clutter 


3.4  Receiver  Noise 

3.4.1  Noise  flindamentals 


Receivers  generate  thermal  noise  which  masks  weak  echoes  being  received  from  the 
radar  transmtesions.  This  noise  is  one  of  the  fundamental  limitatioiw  on  the  radar  range 
and  is  shown  in  the  range  equation  for  a  bi-static  radar,  where  the  transmitter  and  receiver 
are  separated  (see  Appendix  D). 


Most  of  this  noise  originates  in  the  input  stages  of  the  receiver.  The  reason  is  not  that 
these  stages  are  inherently  more  noisy  than  others  but  when  amplified  by  the  receiver’s  frill 
gain,  noise  generated  there  swamps  out  the  noise  generated  further  along.  Since  the  noise 
and  the  received  signals  are  thus  amplified  equally,  in  computing  signal-to-noise  ratios,  the 
factor  of  receiver  gain  can  be  eliminated  by  determining  the  signal  strength  at  the  input 
to  the  receiver  and  dividing  the  noise  output  of  the  receiver  by  the  receiver  gain.  Thus, 
the  receiver  noise  is  commonly  defined  as  noise  per  unit  of  receiver  gain. 


Noise  at  output  of  receiver 

Receiver  noise  = - — - : - - 

Receiver  Gam 


(22) 
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More  commonly,  the  noise  performance  Figure  of  merit  is  known  as  the  ‘noise  figure’, 
Fn-  This  is  the  ratio  of  the  noise  output  of  the  actual  receiver  to  the  noise  output  of  a 
hypothetical,  ‘ideal’  minimum-noise  receiver  providing  equal  gain. 

^  _  Noise  output  of  actual  receiver 
”  Noise  output  of  ideal  receiver  ^ 

Thermal  noise  is  spread  more  or  less  uniformly  across  the  entire  spectrum,  [8].  So,  the 
amount  of  noise  appearing  in  the  output  of  the  ideal  receiver  is  proportional  to  receiver 
bandwidth.  The  mean  power  per  unit  of  receiver  gain,  of  the  noise  in  the  output  of  the 
the  hypothetical  ideal  receiver  is  thus: 

Mean  noise  power  (ideal  receiver)  =  kToB  {W)  (24) 

where  B  is  the  receiver  bandwidth  and  by  convention.  To  is  taken  to  be  290°if ,  which  is 
close  to  room  temperature  and  makes  fcTo  a  round  number  (4  x  10“^^  Ws). 

When  the  internally  generated  noise  is  considerably  greater  than  the  external  noise, 
the  noise  figure,  Fn,  multiplied  by  the  previous  expression  for  mean  noise  power  per  unit 
of  gain  for  an  ideal  receiver  is  commonly  used  to  represent  the  level  of  background  noise 
against  which  target  echoes  must  be  detected. 

Mean  noise  power  (actual  receiver)  =  FnkToB  {W)  (25) 

Typically  though,  the  noise  temperature  of  the  radar  receiver  is  not  the  most  important 
characteristic,  since  choosing  a  low  noise  amplifier  will  sacrifice  other  important  perfor¬ 
mance  characteristics  such  as  dynamic  range,  instantaneous  bandwidth,  phase  amplitude 
stability  and  cooling  in  the  receiver. 

3.4.2  Noise  Statistics 

It  is  important  that  the  noise  modelled  in  this  section  accurately  produces  the  correct 
statistics,  otherwise  the  target  detection  schemes  in  section  4  will  not  work  properly.  To 
this  end,  the  following  hypothesis  for  the  output  signal  is  defined  and  a  derivation  of  the 
Probability  Distribution  Functions  (PDF)  for  each  case  follows: 

Hq  :  noise  only  is  present  ^  n{t) 

Hi  :  signal  and  noise  are  present  ^  s{t)  +  n{t) 

The  additive  zero  mean  white  noise  n{t),  is  spatially  incoherent  and  uncorrelated  with  vari¬ 
ance  cr^.  It  is  added  to  both  magnitude  and  phase  and  hence  has  a  complex  representation. 

Consequently,  the  front  and  rear  signals  can  be  written  as  a  combination  of  real  (in- 
phase)  and  imaginary  (quadrature)  parts,  separated  by  a  phase  of  90‘’. 

^{t)  ~  cos(27r/i„t)  4-  UQ{t)  sin(27r/jnt) 

=  r(t)  cos{2tt fint  - 

or  simply 

=  r(t)  cos(0(t)) 
i^Q{t)  =  r{t)sm{(f){t)) 
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where  fin  =  /c  +  /r  and  fc  +  /f  for  the  rear  and  front  signals  respectively,  r{t)  is  the 
amplitude  of  u{t)  and  the  phase,  <j){t)  =  arctan(^). 


If  we  run  a  test  scenario  where  there  is  no  target,  then  the  in-phase  and  quadrature 
components  are: 

F/{i)  =  ni{t) 

J'Qit)  =  WQ(f) 

and  likewise,  if  we  know  there  is  a  target  present  with  an  amplitude  A: 

2//(f)  =  A  +  ni{t)  nj{t)  —  r{t)  cos{(l>{t))  —  A 
i^Q{t)  =  ngit)  =  r{t)sm{4>{t)). 

where  the  noise  in-phase  and  quadrature  components  nj(t)  and  nQ{t)  are  uncorrelated 
zero  mean  Gaussian  with  equal  variances  a^.  The  joint  PDF  of  the  two  random  variables 
nj-nq  is 

1  /  (r  cos^  —  ^)^  +  (rsin^)^ 

—  exp  I  — ^ - 


27ro’2 


2o'2 


Now,  the  PDF’s  of  the  random  variables  r{t)  and  (p{t),  represent  the  modulus  and  phase 
of  ^{t).  The  Joint  PDF  of  these  two  random  variables  is  given  by 


^  fini;nQ)\J\ 


(27) 


where  the  Jacobian  is: 

J  = 


dni 

dtij 

dr 

dip 

dtiQ 

.  dr 

d<p  . 

L 

COS  (j>  — r  sin  (f 
sin  (f  r  cos  ^ 


with  determinant  |  J|  =  r.  Substituting  this  result  and  equation  26  into  equation  27  gives: 

/(g^) 


r  /  (rcos  (b  —  A)^  -f  (r  sin  Sy 
exp  I  — - 


2-Ka^ 

r 


2Tra^  \  2a^  J  \ 

Then  to  obtain  the  PDF  for  r  alone,  we  need  to  integrate  over 


+  A^ 


2cr2 


exp 


rA  cos  <j> 


r2TT 

fir)  =  /  /(r;#)# 

Jo 


^  exp  - 


+  A‘^\  1  f 


2w 


exp 


^  \  20-2  J  2n  Jq  '  \ 


{rA  cos  <j>) 


where  the  integral  inside  equation  29  is  known  as  the  modified  Bessel  function  of  order 


zero; 


1 
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Thus, 


fir)  =  ^/o 


+  A^'\ 
2^2  J 


(31) 


which  is  the  Rice  PDF.  If  the  scenario  contains  noise  only,  then  ^  =  0  and  equation  31 
becomes  a  Rayleigh  PDF; 

/(r)  =  ^exp(-^)  (32) 

and  if  there  is  a  target  with  the  noise,  then  ^  is  very  large  and  equation  31  becomes  a 
Gaussian  PDF,  with  mean  A  and  variance  cr^: 

The  PDF  for  the  random  variable  (f)  is  found  likewise,  but  not  described  here.  If  we  look  at 
the  noise  statistics  of  the  simulation  over  200  runs,  the  amplitude  distributions  of  target 
with  noise  and  noise  alone  match  the  theoretical  distributions.  One  is  a  Gaussian,  while 
the  other  is  a  Rayleigh  distribution. 


Comparison  of  Probability  Distributions 


Figure  15:  Simulation  probability  distributions 
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3.4.3  Simulation  Results 

In  the  simulation,  the  receiver  noise  N,  used  is  additive  zero-mean  coloured  noise  with 
an  amplitude  determined  by  the  user  defined  noise  power  level.  This  signal  is  added  to 
both  the  front  and  rear  receiver  blocks  in  the  Simulink  simulation.  The  equation  for  the 
doppler  output  signal  can  now  be  updated: 

Sd  =  (Sr  +  N)0(Sf  +  N)  (34) 


Scenario’s  1  and  2  are  shown  in  figures  16  and  17.  They  both  clearly  show  the  effect 
of  the  randoiiiness  over  the  entire  spectrum. 


receiver_nolse_case1 


160 

140 

120 

100 

80 

60 

40 


L . 1 . 1 . 1 . 

. ;  1 . 

. 1 . ::i . j . -t 

_ i _ i _ i - 

_ _ 

_ i _ i _ i - 

80  100  120 
Frequency  (kHz) 


Figure  16:  Scenario  1  with  receiver  noise 


receivef_noise_case2 
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Figure  1 1:  Scenario  2  with  receiver  noise 
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3.5  Complete  Simulation 

The  last  stage  of  the  simulation  is  to  combine  the  three  effects  together.  The  final 
equation  for  the  doppler  output  is: 

Sd  =  (Sr  +  Sm  +  N)  (8>  (Sp  +  Sc  +  N)  (35) 

where  the  different  amplitudes  are  determined  by  the  user.  Simulation  results  are  shown  in 
figures  18  and  19.  The  spectra’s  clearly  show  the  sea-clutter  peak  at  the  correct  frequency, 
the  spread  target  signal  and  random  noise  over  the  entire  spectrum.  When  the  combination 
of  the  sea-clutter  and  multipath  merge  together,  only  the  tip  of  the  target  is  clear. 


mlxed_case1 


Figure  18:  Scenario  1  with  all  three  effects  combined 


mtxed_case2 


Figure  19:  Scenario  2  with  all  three  effects  combined 
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4  Doppler  Detection 

Target  detection  or  ‘lock  on’  is  a  post-launch  ftinction  since  many  missiles  are  not  in  a 
position  to  view  the  target  before  launch.  Those  that  are  able  to  lock  onto  a  target  may 
break  lock  due  to  a  launch  shock,  plume  effects,  or  extremely  high  feedthrough,  which  is 
an  unwanted  DC  component  in  the  spectrum.  Thus,  the  seeker  must  accomplish  target 
detection  at  some  prescribed  time  during  flight. 


Commitment  to  fire  on  a  target  implies  that  a  missile  begiiK  its  task  with  the  knowledge 
that  a  target  exists.  For  a  successfiil  firing,  the  probability  of  finding  a  target  must  be 
95  percent  or  higher.  To  achieve  this,  the  initial  detection  threshold  is  set  relatively  low, 
ensuring  a  high  probability  of  seeing  the  target.  The  probability  of  false  alarm  is  then 
high,  and  a  suitable  detection  scheme  must  be  employed  to  compensate.  After  a  brief 
summary  of  detection  theory,  two  different  doppler  detection  schemes  are  outlined.  The 
first  is  an  optimal  scheme  for  calculating  the  doppler  frequency  and  amplitude,  while  the 
second  is  a  system  employed  in  many  semi-active  systems.  A  simulation  of  the  second 
scheme  has  been  done  and  the  results  are  summarised  in  the  final  section.  In  a  semi-active 
system,  the  target’s  AOA  is  calculated  using  monopulse  principles,  but  that  is  beyond  the 
scope  of  this  report. 


4.1  Detection  Theory 

A  target  is  detected  when  its  amplitude  exceeds  a  threshold  value  r.  The  aim  of  a 
good  detector  is  to  maximise  the  probability  of  detection  for  a  given  probability  of  false 
alarm.  Using  the  hypothesis  and  pdf’s  derived  from  3.4.2,  we  can  define  the  detection  and 
false  alarm  probabilities  for  a  point  X  on  the  spectrum,  as: 

P.  =  P(X>x|H0  =  (36) 

=  (■^)  • 

where  the  variance,  is  the  mean  noise  power  and  the  threshold,  r  can  be  now  calculated 
based  on  a  given  false  alarm  probability, 

T  =  ln(F/o).  (38) 

The  case  where  the  noise  causes  the  target’s  amplitude  to  be  smaller  than  the  threshold 
is  called  a  miss. 

=  P(X  <rm)  =  £  ^  exp  (-££)  ar  (39) 

The  following  figure  shows  the  regions  where  each  of  these  definitions  are  defined. 


Pfa  =  P{X  >  t\Ho) 


Jt 
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MAGNITUDE  {VOLTS/ 0-) 

Figure  20:  Detection  regions 

4.2  Optimal  Detection 

To  find  an  optimal  doppler  detector,  one  approach  is  to  meiximise  a  likelihood  function 
containing  an  accurate  model  of  the  signal.  This  approach  has  been  looked  at,  but  not 
completed.  The  derivation  is  included  in  Appendix  C. 


4.3  Non-Optimal  Detection 

Many  different  techniques  for  doppler  extraction  and  tracking  have  evolved  in  the 
last  30  years.  The  most  common  systems  use  combinations  of  the  following  three  stages: 
narrow-banding,  doppler  frequency  extraction  and  doppler  tracking. 

Narrow- banding 

•  ‘Late  narrow-banding’  segregates  the  wanted  target  signal  in  doppler  frequency  from 
all  other  signals  by  a  narrow  bandpass  filter  at  a  late  stage  in  the  receiver  chain 
incorporated  within  the  main  IF  amplifier. 

•  ‘Early  narrow-banding’  does  the  segregation  as  early  as  possible  in  the  receiver  chain, 
typically  before  the  IF  amplifier. 

Doppler  frequency 

•  ‘Explicit’  doppler  frequencies  are  translated  to  baseband  and  each  doppler  frequency 
is  explicitly  represented  by  an  AC  frequency. 

•  ‘Implicit’  doppler  frequency  remains  always  as  a  shift  on  the  carrier  frequency  and 
its  presence  is  implied  by  the  change  of  frequency. 
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Doppler  tracking 

•  ‘Frequency  locked  loop’  -  The  doppler  frequency  of  the  wanted  target  signal  is 
changed  to  a  second  intermediate  frequency  which  is  compared  with  that  set  by 
a  fixed  frequency  discriminator.  Any  maladjustment  in  frequency  is  sensed  and  used 
to  retune  the  associated  LO. 

•  ‘Phase  locked  loop’  -  A  phase  discriminator  is  used  which  senses  phase  difference. 
This  difference  is  used  to  retune  the  LO  via  an  integrator. 

Target  acquisition  is  accomplished  by  sweeping  the  frequency  of  the  speedgate  LO  over 
the  designated  portion  of  the  doppler  bandwidth.  The  speedgate  LO  is  typically  some  form 
of  voltage-controlled  oscillator  (VCO).  The  linearity  of  its  frequency- voltage  characteristic 
must  be  carefully  controlled  to  allow  precise  sweep  positioning  and  a  constant  sweep  rate 
over  the  full  doppler  frequency  range  of  interest.  Also,  since  the  VCO  forms  part  of  the 
doppler  tracking  circuit,  control  linearity  is  required  to  maintain  a  constant  loop  gain  over 
its  total  frequency  coverage  range.  The  sweep  rate  is  based  on  the  doppler  filter  bandwidth 
and  must  be  slow  enough  to  ensure  signal  buildup  in  the  filter.  The  extent  of  the  sweep 
will  differ  from  system  to  system,  depending  on  the  accuracy  of  the  designation  and  on 
operating  frequency.  Together,  the  rate  and  extent  of  the  sweep  will  determine  the  num¬ 
ber  of  looks  at  the  target  during  the  available  acquisition  time  and  hence  the  cumulative 
detection  probability. 

The  actual  detection  process  consists  of  programming  the  speedgate  LO  with  a  saw¬ 
tooth  of  triangular  sweep  voltage.  When  the  difference  between  the  LO  and  the  target 
doppler  frequencies  equals  the  speedgate  filter  center  frequency,  an  output  is  produced  at 
the  speedgate  discriminator  or  at  a  separate  amplitude  detector.  The  output  signal  is  then 
detected  and  if  it  exceeds  the  detection  threshold,  the  search  is  stopped  for  a  few  tens  of 
milliseconds  while  the  signal  is  examined  to  verify  that  it  is  a  coherent  target  and  not  a 
false  alarm  due  to  noise.  This  process  of  verification  examines  the  signal  in  the  gate  for 
persistence,  since  a  target  will  remain  above  the  verification  threshold  while  noise  will  not. 
Typically,  the  false  alarm  rate  and  the  verification  time  required  for  each  false  alarm  can 
be  optimised  for  a  given  system. 


Once  the  verification  threshold  has  been  passed,  the  speedgate  tracking  loop  is  closed 
and  the  speedgate  is  said  to  be  ‘locked  on’.  A  valid  target  is  then  tracked  in  frequency  using 
the  discriminator,  and  guidance  commands  can  be  extracted  from  it.  A  block  diagram 
representing  the  system  is  shown  in  Figure  21,  [6]. 
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Figure  21:  Alternative  methods  of  target  detection  and  verification 


4.4  Detection  Simulation 


A  simulation  based  on  Figure  21  was  created  in  MATLAB  to  test  the  effect  of  noise, 
multipath  and  sea-clutter  on  doppler  detection.  The  tracking  discriminator  is  not  present 
in  this  simulation,  since  we  are  only  interested  in  the  doppler  detection.  Consequently, 
the  system  is  simplified  and  can  be  represented  by  the  block  diagram  in  Figure  22. 


Doppler 

Spectrum 


Detection 


Figure  22:  Simplified  block  diagram  for  detection  and  verification 

Since  the  worst  case  centre  frequency  of  the  MLC  is  at  the  velocity  of  the  missile,  the 
‘doppler  spectrum’  or  ‘passband’  can  be  determined  as  a  portion  of  the  total  spectrum 
coming  from  the  mixer,  with  the  start  point  set  by  moving  a  fixed  distance  away  from  the 
MLC  and  the  end  point  set  at  a  harmonic  of  the  MLC  frequency.  In  this  simulation,  the 
third  harmonic  was  chosen  as  the  end  point. 
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There  are  a  number  of  important  parameters  in  this  simulation.  These  include  the 
sweep  rate  or  number  of  steps  through  the  filter  bank  over  a  given  time,  the  bandwidth 
of  the  filter,  whether  to  overlap  the  spectrum  from  one  filter  block  to  the  next,  the  detec¬ 
tion  range,  pause  length  and  the  detection  threshold.  The  consequence  of  a  non-optimal 
detection  scheme,  is  that  the  relationship  between  the  false  alarm  probability  and  thresh¬ 
old  calculated  in  section  4.1  no  longer  holds.  The  theoretical  calculation  was  based  on  a 
system  using  the  entire  spectrum  and  not  filtered  sections  of  it.  The  best  analysis  that 
can  be  done  is  a  comparison  of  the  parameters  and  how  an  ideal  simulation  compares  to 
a  realistic  one. 


4.4.1  Doppler  Filter  Bank  and  the  Overlapping  Spectrum 

In  Figure  21,  a  LO  is  used  to  mix  the  doppler  spectrum  down  to  baseband,  so  it  can 
be  used  in  a  single  filter.  For  the  purposes  of  modelling  the  detector,  this  is  simplified 
by  creating  a  new  bandpass  finite  impulse  response  filter  at  each  iteration  through  the 
detection  loop.  The  overall  effect  is  the  same  as  if  a  box  is  drawn  around  the  LO,  mixer 
and  filter. 

In  this  simulation,  the  bandwidth  is  altered  depending  on  the  sweep  rate  and  there  is 
an  option  for  using  a  50%  overlap  between  one  filter  block  and  the  next.  This  increases 
the  sweep  rate  and  minimises  the  loss  of  output  when  a  signal  lies  between  the  center 
frequencies  of  two  filter  blocks.  For  example,  if  there  are  N  filter  blocks  in  a  given  time 
period,  then  there  are  2N  —  1  when  using  a  50%  overlap. 


Figure  23:  Comparison  of  non-overlapping  and  overlapping  filters 


4.4.2  Detection  Range 

The  detection  range  is  required  to  test  ffi,  the  case  where  there  is  a  signal  and  noise 
present.  It  defines  the  bounds  where  a  correct  detection  takes  place  relative  to  the  center 
frequency.  If  a  detection  takes  place  outside  the  detection  range,  it  is  deemed  a  miss. 
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Filtered  Input  Signal 


4.4.3  Pause  Length 

The  pause  length  determines  how  long  the  system  waits  before  deciding  if  a  signal 
contains  a  target  and  not  noise.  If  the  pause  length  is  longer,  there  will  be  less  noise 
that  can  be  mistaken  for  a  target,  but  the  miss  rate  will  be  higher  since  there  is  a  longer 
time  interval  where  a  target  can  not  be  detected.  The  following  figure  has  a  target  with 
noise  present  for  the  first  15  time  periods,  then  noise  only  for  15  time  periods.  There 
are  four  pairs  of  images,  each  pair  having  the  thresholded  only  scenario  on  the  left  and 
the  thresholded  and  delayed  scenario  on  the  right.  If  a  filter  block  is  white,  then  it  has 
exceeded  the  threshold  and  is  a  detection. 

The  target  is  clear  in  each  of  the  scenarios  as  the  constant  line  in  the  bottom  left.  As 
the  pause  length  increases,  more  noise  is  removed  until  only  the  target  is  seen  in  the  final 
scenario  on  the  bottom  right  image. 


4.5  Receiver  Operating  Characteristics 

To  do  a  complete  comparison  of  the  different  parameters,  a  common  technique  is  to 
vary  the  threshold,  and  measure  the  false  alarm  and  detection  rates.  This  is  called  a 
Receiver  Operating  Characteristic  (ROC).  Recall  the  two  hypothesis: 

Hq  :  noise  only  is  present 

Hi  :  signal  and  noise  are  present 

The  axes  of  the  ROC  curve  relate  Hq  to  the  false  alarm  rate  and  Hi  to  the  detection  rate 
using  the  definitions  from  section  4.1. 

The  following  table  describes  the  default  variables  used  in  the  simulation  for  all  pa¬ 
rameters  except  the  one  being  varied. 
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Comparison  of  Delay  Lengths  (Time  vs.  Filter  Block) 
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Figure  25:  Comparison  of  pause  lengths  on  detection 


Table  2:  Detection  parameters 


50%  Overlap  [  Threshold  Pause  (time  periods)  |  Detection  Range  |  Sweep  Rate 
No  I  2  I  3kHz  I  50 


The  following  ROC  curves  were  then  obtained  by  running  both  scenarios  75  times 
for  Ho  and  75  times  for  Hi.  The  sea-clutter  and  multipath  effects  were  included  in  the 
simulation,  while  the  receiver  noise  was  re-simulated  each  time  to  generate  the  randomness 
required. 

The  effect  of  overlapping  filter  banks  caused  an  increase  in  the  false  alarm  rate  for 
both  scenarios.  This  is  a  consequence  of  the  increased  sweep  rate  required  to  compensate 
for  the  overlapping  filter  bank.  The  second  scenario  has  approximately  twice  the  false 
alarm  rate  than  the  first  due  to  the  closer  proximity  of  the  sea-clutter  peak.  Further  work 
could  be  done  to  verify  this  by  changing  the  position  of  the  sea-clutter  and  observing  the 
simulation  results. 
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Ideal  Simulation  ROO 


Multipath  Simulation  ROC 


Figure  26:  ROC  curve  for  varying  overlap  -  scenario  1 


Figure  27:  ROC  curve  for  varying  overlap  -  scenario  2 


The  curves  look  very  similar  in  both  scenarios,  though  the  second  one  has  a  slightly 
higher  false  alarm  rate.  There  does  not  appear  to  be  many  differences  between  the  ideal 
and  the  multipath  cases  indicating  that  this  variable  is  independent  of  the  scenario.  Clearly 
though,  as  the  pause  length  increases,  the  false  alarm  rate  decreases  in  each  case. 
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Figure  30:  ROC  curve  for  varying  detection  range  -  scenario  1 
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Figure  31:  ROC  curve  for  varying  detection  range  -  scenario  2 

All  but  the  multipath  case  in  scenario  1  were  able  to  give  distinctive  results,  as  the 
curves  overlap.  The  trend  is  as  the  sweep  time  increases,  the  false  alarm  rate  gets  higher. 
This  is  due  to  the  bandpass  filter  having  a  smaller  bandwidth  as  the  sweep  time  increases. 
There  is  also  a  slight  increase  in  false  alarm  rate  for  scenario  2.  Prom  these  results  there 
are  a  few  conclusions  that  can  be  drawn: 

•  The  multipath  case  did  not  give  significantly  worse  results,  as  the  broad  sea-clutter 
overlapped  the  multipath  component. 

•  Using  a  50%  overlap  in  the  filter  bank  gives  a  higher  false  alarm  rate. 

•  A  longer  threshold  pause  gives  a  lower  false  alarm  rate  and  is  independent  of  the 
scenario. 

•  The  trend  for  the  detection  range  shows  that  larger  detection  ranges  give  lower  false 
alarm  rates  until  an  upper  limit  is  reached. 

•  A  greater  sweep  time  (filter  bank)  gives  a  higher  false  alarm  rate. 
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Ideal  Simulation  ROC 


l^ultipatH  Simulation  R^>C 


Figure  32:  ROC  curve  for  varying  sweep  rates  -  scenario  1 


Multipath  Simulation  ROC 


Figure  S3:  ROC  curve  for  varying  sweep  rates  -  scenario  2 

5  Conclusions  and  Future  Work 

The  Simulink  simulation  worked  well  and  demonstrated  the  effects  of  multipath,  sea- 
clutter  and  receiver  noise  on  an  ideal  simulation.  It  is  clear  that  mixing  the  multipath  and 
sea-clutter  caused  the  spectrum  to  spread  widely  and  this  in  turn  gave  some  interesting 
detection  results.  There  are  a  number  of  improvements  that  could  build  on  this  section: 

•  Multipath  with  a  reduced  amplitude  could  be  added  to  the  front  signal. 

•  A  small  randomness  in  the  frequency  offset  for  the  multipath  signal  could  be  used  to 
more  accurately  model  real  life.  This  would  add  extra  frequency  components  during 
the  convolution. 

•  A  comparison  on  different  amounts  of  frequency  spread  could  be  included. 

•  The  models  for  multipath  and  sea-clutter  could  be  improved  to  more  accurately 
model  different  sea  states  and  other  atmospheric  conditions. 
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The  optimal  detection  scheme  needs  some  more  work  to  provide  an  alternative  means 
to  determine  the  doppler  frequency  and  amplitude  of  the  filtered  spectrum.  The  non- 
optimal  detection  scheme  gave  good  results  and  led  to  a  number  of  observations.  These 
are  repeated  for  completeness; 

•  The  multipath  case  did  not  give  significantly  worse  results,  as  the  broad  sea-clutter 
overlapped  the  multipath  component. 

•  Using  a  50%  overlap  in  the  filter  bank  gives  a  higher  false  alarm  rate. 

•  A  longer  threshold  pause  gives  a  lower  false  alarm  rate  and  is  independent  of  the 
scenario. 

•  The  trend  for  the  detection  range  shows  that  larger  detection  ranges  give  lower 
false  alarm  rates  until  an  upper  limit  is  reached.  More  extensive  simulations  of  this 
variable  could  give  more  conclusive  results  and  be  used  to  optimise  the  detection 
range  for  this  system. 

•  A  greater  sweep  time  gives  a  higher  false  alarm  rate,  due  to  the  decreasing  bandwidth 
of  the  bandpass  filter.  Again,  further  simulations  could  be  used  to  find  an  optimal 
sweep  time  for  this  system. 


The  detection  simulation  could  be  improved  by  increasing  the  length  of  time  that  each 
hypothesis  was  tested  over  and  increasing  the  number  of  points  on  the  ROC  curves.  Fur¬ 
ther  simulations  could  also  be  used  to  optimise  the  detection  range  and  the  sweep  time 
for  this  system.  There  are  also  a  few  alternative  detection  schemes  that  could  be  used  to 
compare  with  this  one,  such  as  the  ‘n  out  of  m  detector’,  [6],  detectors  based  on  a  ‘constant 
false  alarm  rate’  and  schemes  which  incorporate  both  detection  and  tracking. 


Finally,  putting  the  seeker  in  the  real  world  would  involve  the  target  moving  in  different 
directions  and  at  different  velocities.  The  seeker  would  need  to  compensate  by  changing 
its  trajectory  and  velocity.  Both  of  these  changes  would  vary  the  doppler  measured  at  the 
front  receiver.  In  the  future,  this  model  could  be  introduced  into  a  larger  system  where 
there  are  moving  targets  and  existing  sea-clutter  and  multipath  models.  Testing  the  effect 
under  these  dynamic  conditions  is  a  much  larger  project,  but  could  easily  extend  this  work. 
The  ultimate  outcome  would  then  be  to  verify  ‘Hardware  In  the  Loop’  simulations  with 
the  measured  output  doppler  and  AO  A. 


31 


DSTO-TR-1606 


References 

1.  W.V. Andrew,  The  effects  of  Multi-path  Propagation  on  Low-Altitude  Detect  Pulsed- 
Doppler  RADAR  Systems,  December  1990,  Arizona  State  University. 

2.  W.Lippincott,  Sea  Surface  Multipath  Effects  on  Ship  Radar  Radiated  Power  Determina¬ 
tion,  Advanced  Systems  Technology  Branch,  Space  Systems  Development  Department. 

3.  R.M.Smith,  J.Y.Yee,  C.S.  An  and  A.L.Haun  Simulation  of  Multipath  for  Semi-active 
Missiles,  AGARD  CP-dTS,  Turkey  1990, 

4.  D.Lewis  Principles  of  Naval  Architecture,  1983,  Society  of  Naval  Architects  and  Marine 
Engineers,  August  1997. 

5.  C.  Beard  and  I.Katz  The  Dependence  of  Microwave  Radio  Spectra  on  Ocean  Roughness 
and  Wave  Spectra,  IVansactions  on  Antennas  and  Propagation,  1957 

6.  M.Skolnik,  "‘RADAR  Handbook  -  Second  Edition,  1990,  McGraw-Hill  Inc. 

7.  D.A.James,  ‘Radar  Homing  Guidance  for  Tactical  Missiles,  1986,  Royal  Military  Col¬ 
lege  of  Science  Shrivenham. 

8.  G.W.Stimson,  Introduction  to  Airborne  RADAR,  1983,  Hughes  Aircraft  Co.,  El  Se- 
gundo,  CA. 

9.  R.Smith,  Multipath  Math  Model  and  Implementation  Guide  -  Revision  No.  1,  July  29, 
1986,  Missile  Software  Branch,  China  Lake,  California  USA. 

10.  B.Moran  Detection,  Estimation  and  Classification,  MSIP  Lecture  notes,  March  2000. 

11.  B.Mahafza  Radar  Systems  Analysis  and  Design  Using  MATLAB,  Chapman  and 
Hall/CRC,  2000 


32 


DSTO-TR-1606 


Appendix  A  Ideal  Simulation  Verification 


The  Simulink  model  used  for  the  simulation  is  shown  in  Figure  Al.  The  low  pass  filter 
is  implemented  in  Simulink  by  first  creating  it  in  MATLAB,  then  passing  its  frequency 
response  to  be  multiplied  in  the  frequency  domain. 


Figure  Al:  Simulink  model 


Using  this  model,  scenario  1  was  run  with  a  rear  doppler  of  -30kHz  and  a  front  doppler 
of  50kHz.  With  the  carrier  reference  set  to  IMHz  rather  than  lOGHz  for  the  convenience 
of  demonstration,  the  frequency  of  the  front  signal  should  be  lMHz-l-50kHz— 1.05MHz  and 
the  frequency  for  the  rear  signal  should  be  lMHz-30kHz=0.97MHz.  The  front  and  rear 
signals  are  shown  in  Figure  A2: 
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Figure  AS:  Ideal  scenario  1  repeated 
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Appendix  B  Sea-clutter  Frequency  Calculation 

To  calculate  the  sea-clutter  center  frequency  shift,  a  Gaussian  shaped  spectrum  is  ap>- 
proximated  by  the  Fourier  Transform  of  a  cosine  function,  where  the  mean  is  approximated 
by  the  frequency  of  the  cosine.  The  exact  frequency  was  then  determined  by  the  following 
signal  definitions  in  the  time  domain: 

Spit)  =  sin(27r(/c  + /f)^) 

SR{t)  =  sin(27r(/c  -1-  fR)t) 

Sc{t)  =  sin(27r(/c  +  fci)t) 

where  the  new  front  signal  is  now  represented  as  in  equation  19  and  the  mixing  is  a 
multiplication,  since  the  signals  are  represented  in  the  time  domain.  The  final  spectrum 
becomes: 

Snit)  =  SR{t)[SF{t)  +  Sc{t)] 

=  i  cos(27r(/F  -  fR)t)  +  i  cos(27r(/cz  -  fR)t)  + 

^  cos(27r(2/c  +  fR  +  fF)t)  +  i  cos(27r(2/c  +  /f  +  fcl)t)  (Bl) 

where  the  desired  sea-clutter  frequency  is  now  shifted  by  the  rear  doppler  frequency. 
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Appendix  C  Optimal  Detection 


£,Prom  the  complete  simulation  in  section  3.5,  we  can  define  a  signal  s(/),  in  the 
frequency  domain  with  all  three  effects  added,  sea-clutter,  multipath  and  noise.  Math¬ 
ematically,  s{f)  is  multiplied  by  a  Gaussian  with  unknown  amplitude  ay,  mean  /r  and 
known  variance  <j|.  and  contains  additive  zero  mean  white  noise,  n(/)  with  variance  cr^; 


*(/)  =  orexp  (  -  ^ 


(Cl) 


The  probability  density  of  s(/)  is  given  by 


/  I  if  ~  fTY 

s(/)-OTexp( - 


and  the  likelihood  function  is  obtained  by  taking  the  log  of  the  probability  density: 

/ri  =  E  |»(/)  -  I 

If  we  differentiate  this  with  respect  to  the  two  unknowns,  ut  and  fx,  we  get: 

U-hf' 


-^E 

ITiT^  Z— ✓ 


dax  'Ko 

dLs  _  2 
dfx 


/ 


s(/)  -arexp 


s{f)  -  ax  exp  (  - 


2<t| 


.  U-hf 

"’‘p  I  — s— 


(/  -  hY 

2a% 


ay  exp 


{f-fTf\  if-fT 


2(t|, 


2(t|, 


(05) 


If  these  two  equations  are  then  equated  to  zero,  they  become 

if-h? 


E 

/  ■- 


s{f)  -  ax  exp 


2(Tj. 


E 

/ 


,(/)-<.Texpl - ^ 


(CT) 


Equation  C6  can  then  be  rearranged 


i:»(/)-p  (-^) = -E-p  (-^)  (C8) 


to  show  that  the  RHS  represents  a  Gaussian  probability  density  function,  with  mean  fx 

U-hf 


and  variance  <j|.. 


:  exp 


27rcri 


2(t| 


Hence,  equation  C8  becomes 


^s(/)exp  ^  =  «r^27r<Tf 


(CIO) 
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and  can  be  rearranged  to  give  an  expression  for  ay. 

Now,  if  equation  C7  is  rearranged  in  a  similar  way, 

"  if-fT?' 


(Cll) 


X]  ^(/)(/  -  fT)  exp  (-  —2/^^  )  (/-/r)  (C12) 


/  ^  ^  f 

and  we  use  the  fact  that  E{f  —  fx)  =  0.  Then, 


f 

and  equation  C12  becomes 


(/  -  ItY 

2(t|, 


(/-/r)  =  0 


E''(/)(/-/riexp(-R^)=0. 


(/  -  hy 

2a^ 


The  next  step  is  to  define 

aif  -  It)  =  exp 
and  notice  that  equation  C14  is  the  same  as 

-  /t)  ~  f  s{f)g'{f  -  fT)df  =  0 

Using  integration  by  parts,  this  becomes 

s{f)gU  -  fr)  -  f  s'{f)g{f  -  fT)df  =  0 
/  J  f 

where  the  first  part  tends  to  0.  Thus 

s\f)gif  -  /r)d/  =  0 

Now,  using  Parsevals  Identity: 

I  HP(f)mG{f))df=^  J  mg{t)dt 

we  can  take  the  Fourier  Transform  of  this  expression,  where 

PU’if)}  =  ts{t) 

Half-h)}  = 

r  /t2+2i 

=  e-2^UTt(7exp 

Thus,  equation  C17  becomes 


j  Ht) 


exp 


aH'^ 


g-27r/rt^^  =  0 


(C13) 


(C14) 


(CIS) 

(C16) 

(C17) 

(CIS) 

(C19) 

(C20) 

(C21) 


At  this  stage,  a  numerical  solution  in  MATLAB  is  required  to  solve  this  equation. 
Work  on  this  method  is  still  in  progress. 
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Appendix  D  Range  Equation 

Here  is  the  range  equation  for  a  bi-static  radar,  as  used  in  the  semi-active  case  in  this 
report. 


(RlRl) 


max 


PjcrGjAffi 

{4Tr)'^LSmin 


(Dl) 


where, 

Ri  =  range  illuminator-target  (m) 
i?2  =  range  target-missile  (m) 

Pi  =  ill  uminator  power  (W) 
a  =  radar  cross  section  (m^) 

G]  —  illuminator  gain  (dB) 

Afn  =  effective  aperture  of  the  receiving  antenna  (m^) 
L  =  loss  factor  for  miscellaneous  Losses  {dB) 


and,  Smin  —  Fn{kToB){S  i  iV)min>  where 


{S  •  N)min  — 
Fn  = 
B  - 
k  = 
To  = 


minimum  detectable  signal  level  at  the  receiver  output  (dB) 
receiver  noise  figure 
receiver  bandwidth  (Hz) 

Boltzmann’s  constant  1.38  x  (Ws/^K) 

absolute  temperature  of  the  resistor  representing  the  external  noise  (^K) 
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